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ABSTRACT 

We compute the relic abundance of the right-handed sneutrinos in the supersymmetric 
Fo-term model of hybrid inflation. As well as providing a natural solution to the /i- and 
gravitino overabundance problems, the Fo-term model offers a new viable candidate to 
account for the cold dark matter in the Universe: the lightest right-handed sneutrino. 
In particular, the Fo-term model predicts a new quartic coupling of purely right-handed 
sneutrinos to the Higgs doublets that thermalizes the sneutrinos and makes them annihilate 
sufficiently fast to a level compatible with the current cosmic microwave background data. 
We analyze this scenario in detail and identify favourable regions of the parameter space 
within the framework of minimal supergravity, for which the lightest right-handed sneutrino 
becomes the thermal dark matter, in agreement with WMAP observations of cosmological 
inflation. Constraints derived from direct dark matter searches experiments are presented. 
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1 Introduction 



Hybrid inflation [1], along with its supersymmetric realizations [2-5], remains one of the 
most predictive and potentially testable scenarios of inflation that have been suggested so 
far. Hybrid inflation is predictive and testable, in the sense that the inflaton dynamics 
is mainly governed by a few renormalizable operators which might have observable impli- 
cations for laboratory experiments. In such a scenario, inflation terminates through the 
so-called waterfall mechanism, which is triggered, when the inflaton field <f> passes below 
some critical value <p c . From that point on, another field X, called the waterfall field, held 
fixed at origin initially, quickly rolls down to its true vacuum expectation value (VEV) and 
drastically modifies the slow-roll form of the 0-potential, thereby ending inflation. 

In supersymmetric theories, the required form of the hybrid inflationary potential may 
originate from either the F-terms of the superpotential or from a large Fayet-Iliopoulos (FI) 
D-term [6], usually induced by some anomalous local U(l) symmetry within the context 
of string theories. In both the F- and D-term hybrid inflation, the slow-roll slope of 
the potential may come either from supergravity (SUGRA) corrections [2] and/or from 
radiative effects [3-5]. 

Recently, a new supersymmetric hybrid inflationary model was proposed in [7] and 
studied in detail in [8]. The model realizes F-term hybrid inflation and includes a sub- 
dominant non- anomalous FI D-term that arises from the U(l)x gauge symmetry of the 
waterfall sector. It has therefore been called the Fp-term model of hybrid inflation, or in 
short, the Fp-term model. The F^-term model can naturally accommodate the currently 
favoured red-tilted spectrum with n s — 1 —0.037 [9], along with the actual value of the 
power spectrum of curvature perturbations, P-r ~ 4.86 x 10~ 5 [10], and the required number 
of e-folds, M e ~ 55 [11]. 

The presence of the FI term in the F D -term model is necessary to approximately break 
a F-parity that governs the waterfall sector. The approximate breaking of the F-parity 
gives rise to late decays of the superheavy waterfall-sector particles that are produced just 
after inflation during the preheating epoch [12,13]. These waterfall particles have masses 
of the Grand Unified Theory (GUT) scale and can dominate the energy density of the 
Universe, provided the inflaton coupling k to the waterfall sector is not too suppressed, 
i.e. for values of k > 10~ 3 . Then, the late decays of the GUT-scale waterfall particles 
produce an enormous entropy that can reduce the gravitino abundance Yq well below the 
limits imposed by big bang nucleosynthesis (BBN), i.e. Yq < 10~ 15 [14]. In this way, the F D - 
term model provides a viable solution to the gravitino overabundance problem [8], without 
the need to unnaturally suppress all renormalizable inflaton couplings to the particles of 
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the Minimal Supersymmetric Standard Model (MSSM) sector, below the 10 6 level. 

Another interesting feature of the Fo-term model is that the p-parameter of the 
MSSM can be generated effectively by the superpotential operator XSHuH^, when the 
scalar component of the inflaton chiral multiplet S receives a non-zero VEV after the spon- 
taneous symmetry breaking (SSB) of the local U(l)x symmetry of the waterfall sector [15]. 
Moreover, the inflaton superfield S couples to the right-handed neutrino superfields Aq^, 
via the superpotential coupling ^pijSNiNj, with i, j = 1, 2, 3. Hence, the inflaton VEV will 
produce an effective Majorana mass matrix as well [7,16]. As a consequence, the resulting 
heavy Majorana neutrinos are expected to have masses of order p. If pij is approximately 
SO (3) symmetric, i.e. p^- ~ p 1 3 , a possible explanation of the observed baryon asymmetry 
in the Universe (BAU) may be obtained by thermal electroweak-scale resonant leptogenesis, 
in a way independent of any pre-existing lepton- or baryon-number asymmetry [17]. 

Even though the Fp-term model violates explicitly the lepton number (L) by AL = 2 
superpotential operators, it conserves F-parity. Hence, the lightest supersymmetric particle 
will be stable and so will potentially qualify as a candidate for the cold dark matter (DM) in 
the Universe. Most interestingly, the F^-term model provides a new candidate for the cold 
DM. This is the lightest right-handed sneutrino (LRHS), which may possess thermal relic 
abundance [8] for relatively large values of the aforementioned superpotential couplings A 
and p, i.e. for A, p > 1CT 2 . This should be contrasted with what is happening in standard 
seesaw extensions of the MSSM, where Aq )2 ,3 have only bare Majorana masses. Because 
the small neutrino Yukawa couplings are the only possible interactions of sneutrinos with 
matter in these models, purely right-handed sneutrinos turn out to be non-thermal and 
tend to overdose the Universe by many orders of magnitude [18,19]. It is therefore difficult 
for the LRHS to be a thermal DM in seesaw extensions of the MSSM with bare Majorana 
masses. 

In this paper we analyze in detail the relic abundance of the right-handed sneutrinos 
in the supersymmetric F^-term model of hybrid inflation. In this model, the F-term of the 
inflaton superfield, Fs, gives rise to the new quartic coupling, ^Ap N*N*H u H d , in the scalar 
potential, which involves the right-handed sneutrinos Ni j2 ,3 and the Higgs doublets H u d . 
As mentioned above, unless the couplings A and p are too small, the new quartic coupling 
will be sufficiently strong to thermalize the sneutrinos and make them annihilate to a level 
compatible with the current cosmic microwave background (CMB) data [9], from which 
the DM component of the Universe was found to be 

tt DM h 2 = 0.1099 ± 0.0062. (1.1) 

The central goal of our analysis is to delineate the parameter space within the context of 
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minimal supergravity (mSUGRA), for which the LRHS is the thermal DM. In addition, we 
consider the constraints obtained by WMAP observations related to cosmological inflation. 
Finally, we present numerical estimates of the scattering cross-section of the LRHS with 
nuclei that will be relevant to direct DM searches in present and future experiments. 

After this introduction, the paper is organized as follows: in Section 2 we present 
the basic structure of the F^-term model and briefly review the solution to the gravitino 
overabundance problem. Moreover, in the same section we derive the constraints imposed 
on the theoretical parameters by cosmological inflation. In Section 3 we perform a detailed 
study of the relic abundance of the LRHS and offer numerical estimates of representative 
scenarios within the mSUGRA framework. We also present numerical estimates for the 
scattering cross-section of the LRHS with the nucleon, indicating the presently achieved 
and future sensitivity of the current and projected experiments for DM searches, such as 
CDM-II, SuperCDMS and XenonlT. Finally, we summarize our conclusions in Section 4. 



2 The Fo-Term Model of Hybrid Inflation 

In this section we first outline the basic structure of the F^-term model of hybrid inflation. 
Then we briefly review how the gravitino abundance can be solved within the Fo-term 
model. Finally, we present the constraints on the theoretical parameters that are imposed 
by CMB data pertinent to inflation. A more detailed discussion of all the above issues may 
be found in [8]. 

2.1 The Model 

The Fo-term model may be defined through the superpotential 
W = KSfx^ - M 2 ) + \SH u H d + ^fSN i N j + h^L^Nj + , (2.1) 

where S is the gauge-singlet inflaton superfield and Xi >2 is a chiral multiplet pair of 
the so-called waterfall fields which have opposite charges under the U(l)x gauge group, 
i.e. Q(Xi) = —Q(X 2 ) = 1. In addition, indicates the MSSM superpotential without 

the /i-term, 

^m7s°m = KQiHjjj + h d ij H d Q i D j + htHMEt . (2.2) 
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Within the SUGRA framework, the sector of soft supersymmetry (SUSY) breaking (SSB) 
derived from (2.1) is given by 

- £ soft = MlS*S+M 2 ^N*N i + (kA k SX 1 X 2 +\A x SH u H d +^ A p SN i N i -nasM 2 S + H.c. ) , 

(2.3) 

where M§, M^, A K ^ P and as are soft SUSY-breaking mass parameters that are all typically 
of order MgusY ~ 1 TeV. In addition, the F^-term model contains a FI D-term, — ^gm^D, 
associated with the U(l)x gauge symmetry of the waterfall sector. The latter gives rise to 
the D-term potential 

V D = j (lYxl 2 - \X 2 \ 2 - m 2 Fi y , (2.4) 

where g is the U(l)x gauge-coupling constant. The FI mass parameter m-pi is subdominant 
with respect to the superpotential tadpole mass M, i.e. m F j/M < 1CT 5 . 

An interesting feature of the F^-term model is the generation of an effective //-term 
of the required order Msusy after the SSB of U(l)x- To see this, let us neglect the VEVs 
of H U:d next to the large VEVs of the waterfall fields Xi >2 : (Xi j2 ) — M. To a good 
approximation, the VEV of S may then be determined by the following part of the potential: 

V s = \F Xl \ 2 + \F X2 \ 2 + M 2 S S*S + [KM 2 (A K -as)S + H.c] , (2.5) 

where we have set the waterfall fields X 12 to their actual VEVs. Substituting the F-terms 
of the waterfall fields, 

F Xh2 = kS(X 2W ) = kMS, (2.6) 

into (2.5), we obtain 

V s = (2k 2 M 2 + M$)S*S + KM 2 (A K -a s )S + H.c. . (2.7) 

It is then not difficult to derive from (2.7) that at the present epoch of the Universe, the 
inflaton field, S, acquires the non-zero VEV 

(S) = y k \A k - a s \ + 0(M 2 VSY /M) , (2.8) 
in the phase convention that (S) is positive. Equation (2.8) implies the effective /z-term 

H = \(S) *i ^\A K -a s \ . (2.9) 

If A ~ k, the size of //-parameter is of order Msusy, as required for a successful electroweak 
Higgs mechanism. 

In addition to the generation of an effective /i-parameter, the third term in (2.1), 
| Pij S NiNj, gives rise to an effective lepton-number-violating Majorana mass matrix, 
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i.e. Mat = PijVs- If we assume that p^ is approximately SO(3) symmetric, i.e. p^ ~ p I3, 
one obtains 3 nearly degenerate right-handed neutrinos Ai i2 ,3, with mass 



m N = pv s ■ 



(2.10) 



If the couplings A and p are comparable, then the /z-parameter will set the scale for the 
SO(3)-symmetric Majorana mass m^, i.e. ~ p [7]. Evidently, this will lead to a scenario 
where the singlet neutrinos iV li2 ,3 have TeV or electroweak-scale masses. This opens up 
the possibility of directly detecting these singlet Majorana neutrinos through their lepton- 
number violating signatures at the LHC [20] or ILC [21]. Furthermore, in the Fp-term 
model the BAU could be explained by thermal electroweak-scale resonant leptogenesis [17]. 

2.2 Solution to the Gravitino Overabundance Problem 

The FI mass term mpi plays a key role in providing a viable solution to the gravitino 
overabundance problem in the F^-term model, without the need to unnaturally suppress 
all the inflaton couplings k, A and p below the 10~ 6 level [7,8]. 

In detail, the presence of m-pi explicitly breaks an unwanted discrete symmetry that 
arises from the permutation of the waterfall fields: X± <-> X 2 . If m-pi was absent, the 
permutation symmetry would remain exact even after the SSB of the U(l)x- This would 
act like parity and was therefore termed D-parity in [7]. As a consequence of D-parity 
conservation, the D-odd waterfall particles of mass gM would have been stable, and if 
abundantly produced during the preheating epoch [12,13], they could overdose the Universe 
at late times. 

To avoid this undesirable situation, we introduce a small but non-zero FI term mpi- In 
this case, the D-odd waterfall particles will have forbidden decays to two Z)-even inflaton- 
related fields of mass kM, induced by the FI term. To kinematically allow for such decays, 
we assume that k < g/2, where g is the value of the U(l)x coupling constant at the GUT 
scale. The late decays of the D-odd waterfall fields will then reheat again the Universe 
at temperature T g , and so release enormous entropy that might be sufficient to reduce the 
gravitino abundance Yq below the BBN limits. More explicitly, after the Universe passes 
through a second reheating phase, the gravitino abundance may be estimated by [8]: 



Hence, for second reheat temperatures T g ~ 1 TeV and inflaton couplings k > 10 2 , the 
strict constraint Yq < 10~ 15 , for < 500 GeV, can be comfortably met. 
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To determine the second reheat temperature T g , we may use the standard freeze-out 
condition T g = H(T g ), where 

- 4 § < 2 - 12 > 

is the decay rate of the D-odd particles and 

is the Hubble expansion parameter in the radiation dominated era of the Universe and 
mpi = 2.4 x 10 18 GeV is the reduced Planck mass. In particular, for a fixed given value 
of T g , we may infer the required size of the FI mass term m^i [8]: 



05\ 3/4 / T g \ 1/2 AO 16 GeV \ 1/4 



8.4X10-X - 77^777 I 77 ■ (2-14) 



— K5 

As can be seen from (2.14), for T 9 ~ 1 TeV, it should be m-pi/M ~ 10~ 6 , so the FI mass 
term tofi needs be much smaller than M. Detailed discussion of how such an hierarchy 
can be naturally achieved within the SUGRA framework may be found in [8]. 



2.3 Constraints from Cosmological Inflation 

Here we recall the constraints derived in [8] on the Fo-term model from cosmological 
inflation. In fact, there are three constraints that need to be considered. 

The first constraint arises from the requirement of solving the horizon and flatness 
problems of the standard Big-Bang Cosmology. According to the inflationary paradigm, 
these problems may naturally be solved, if our observable Universe had an accelerated 
expansion of a number of 50-60 e-folds. In the slow-roll approximation, the number of 
e- folds, Af e , may be calculated by [11] 

K = — / #^ - 55, (2.15) 

m Pl J4>cnd V inf 

where = \[2 Re S is the inflaton field and V^f is the Fp-term inflaton potential that can 
be found in Section 2.1 of [8]. We will always denote differentiation with respect to <fi with 
a prime on V\ n {. Moreover, ex i t is the value of 0, when our present horizon scale exited 
inflation's horizon, whilst en( j is its value at the end of inflation. Specifically, the field 
value en d may be determined from the condition: 

max{e(0 cnd ), |r7(0 cnd )|} = 1, (2.16) 
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with 



m 



n 1 ^\ n - m 2 Ylk (9 17) 



The other two inflationary constraints come from the so-called power spectrum P-ji of 
curvature perturbations and the spectral index n s . The square root of the power spectrum, 
, is given by 

pl/2 1 Vinf 2 (0cxit) (9 1%) 

n ' 2^3^ l^ f (0 cxit )| ' 1 J 

This prediction must be compared with the result obtained by a 3-years WMAP analysis 
of CMB data [10], 

P% 2 ~ 4.86 x 10" 5 . (2.19) 
Moreover, in the slow- roll approximation, the spectral index n s is given by [11] 

n s = l-6e(0 cxit ) + 2r](<f) cxit ) ~ 1 + 2??(0 cxit ), (2.20) 

where the parameter e is negligible in the Fj)-term model. Recently, after analysing its 
data collected in the last 5 years, WMAP has reported the value for the spectral index [9]: 



1 = -0.0371*;™ . (2.21) 



This result slightly favours a red-tilted spectrum and is consistent with scale invariance at 
the 2.64 a confidence level. 

Given the three constraints (2.15), (2.19) and (2.21), and assuming that all inflaton 
couplings are equal, i.e. k = A = p, one obtains within mSUGRA the upper bound [8] 

k < 2x 10~ 2 . (2.22) 

On the other hand, the inflationary scale M is close to the GUT scale, i.e. M ~ 10 16 GeV, 
when k reaches its upper bound imposed by inflation. For an inflaton sector that realizes 
a next-to-minimal Kahler potential with a negative Hubble-induced mass term for S [22], 
the upper limit on k may be slightly relaxed to [8] 

k < 3.2 x 10~ 2 , (2.23) 

whilst M decreases to M ~ 0.5 x 10 16 GeV. 

It is important to properly translate the upper bounds (2.22) and (2.23) on k obtained 
at the inflationary scale M into the respective ones on A and p for the soft SUSY-breaking 
scale Msusy- As we will see more explicitly in the next section, it is the product Ap 
evaluated at the scale M S usy that controls the strength of annihilation of the LRHSs into 
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the Higgs fields and other SM particles. Even though the renormalization group (RG) 
evolution of p from M to Msusy may be ignored, as p(M) pa p(Msusy), this is not the case 
for the coupling A. Neglecting gauge and small Yukawa couplings of order 10" 1 , the RG 
equation for A is given by [23] 

= + 1*0 • (2 ' 24) 

where t = \n(Q 2 /M S vsy)- Assuming that the RG evolution is dominated by the top-quark 
Yukawa coupling h t) the solution to (2.24) is easily found to be 

CM \ 3/l ?/(167r 2 ) 

A(Msusy) = A(MW— pa 0.57 xA(M). (2.25) 

To obtain the last result in (2.25), we assumed that h t ~ 1 and M S usy/^ ~ 10~ 13 . 
Then, starting with the boundary condition A = k at the inflationary scale M, the RG 
running (2.25) of A implies the upper limits: 

A(Msusy) £ 1.14 x 10" 2 , A(M SUSY ) £ 1-82 x 10~ 2 , (2.26) 

for an inflaton sector with a minimal and a next-to-minimal Kahler potential, respectively. 

In addition to constraints from inflation, one may also get constraints on the size 
of M from cosmic strings that arise due to the SSB of the local U(l)x symmetry. For 
values of k ~ 10~ 2 of our interest, this implies that one must have [8] M < 0.5 x 10 16 GeV. 
This constraint may be a bit restrictive for the mSUGRA model, but it can be completely 
avoided if the waterfall sector realizes an SU(2)x gauge symmetry instead of U(l)x, whose 
SSB generates no topological defects [8]. Consequently, we will conservatively consider the 
limits stated in (2.22), (2.23) and (2.26) when implementing inflationary constraints on the 
relic abundance of the LRHS in the next section. 



3 Right- Handed Sneutrino as Thermal Dark Matter 

In the F^-term hybrid model -R-parity is conserved, even though the lepton number L, 
as well as B — L, are explicitly broken by the Majorana operator ^pSNN. We note 
that all superpotential couplings either conserve the B — L number or break it by even 
number of units. Since -R-parity of each superpotential operator is determined to be R = 
(_1)3(b-l) _ _|_^ p D _^ erm hybrid model conserves .R-parity. As a consequence, the LSP 
of the spectrum is stable and can be a viable candidate for Cold Dark Matter (CDM). As an 
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extension of the MSSM, our model can accommodate the standard SUSY CDM candidates, 
such as the lightest neutralino. Because of the connection between the Higgs and neutrino 
sectors, on the one hand, and inflation, on the other, it is very interesting to explore the 
possibility of having a right-handed sneutrino as LSP in order to solve the CDM problem. 
As we will see in Section 3.3, this renders the F^-term model much more constrained, 
leading to sharp predictions for scattering cross-sections relevant to experiments of direct 
searches for CDM. 




3.1 Sneutrino Mass Spectrum 

Before calculating the sneutrino relic abundance in our model, we first observe that light 
right-handed sneutrinos may easily appear in the spectrum. Ignoring the terms proportional 
to the small neutrino- Yukawa couplings, the 6x6 right-handed sneutrino mass matrix M 2 - 
is given in the weak basis (-/Vi,2,3, N*2 3 ) by 

p 2 v\ + M 2 ~ pA p v s + p\v u v d \ 
P A* p v s + p\v u v d p 2 v 2 + M 2 ~ ) ' 

where vg = (S), v u ^ = (H u ^). Moreover, M 2 ~ is the soft SUSY-breaking mass matrix 
associated with the sneutrino fields and A p is the sneutrino trilinear coupling matrix. In 
general, A4 2 ~ is diagonalized by a unitary matrix such that 

4 M % U n = dia § ( m \ ."»#„.■■■. m % 6 ) > ( 3 -2) 
where the sneutrino masses are ordered, such that < m^ 2 < ... < . Neglecting 
the possible flavor structure contained in the 3x3 matrices M~ and A p , the sneutrino 
spectrum will then consist of 3 light (heavy) right-handed sneutrinos with masses 

All mass terms in (3.3) are (9(100-1000) GeV, so a proper choice of model parameters can 
accommodate a LRHS to act as LSP. Unless the trilinear coupling A p is small compared 
to /i, the off-diagonal elements in (3.1) will induce a sizeable mixing between the heavy 
and light right-handed sneutrino states, suppressing the light masses to values smaller than 
(/i 2 + M^) 1 / 2 . This will be demonstrated in our discussion of the numerical results in 
Section 3.3, where the Fo-term model is embedded within the mSUGRA framework. 



3.2 Sneutrino Annihilation and Relic Density 

Right-handed sneutrinos as CDM were considered in [18] in the context of the MSSM with 
right-handed neutrino superfields Ni and bare Majorana masses MjjNiNj. This analysis 



10 



shows that thermal right-handed sneutrinos have rather high relic abundances and will 
generally overdose the Universe. The reason is that because of the small neutrino Yukawa 
couplings h"j, the self- and co-annihilation interactions of the sneutrino LSP with itself and 
other MSSM particles are rather weak. These weak processes do not allow the sneutrino 
LSP to stay long enough in thermal equilibrium before its freeze-out temperature, such 
that its number density gets reduced to the observed value f2oM^ 2 ~ 0.11 [9] [cf. (1.1)]. 
In fact, the predicted values for £~2dm^ 2 turn out to be many orders of magnitude larger 
than I. Instead, right-handed sneutrinos can be viable thermal DM candidates in the 
MSSM if they significantly mix with left-handed sneutrinos, either by increasing the SUSY- 
breaking trilinear couplings [24] 1 , or by lowering the right-handed neutrino mass scale [25]. 
Alternatively, right-handed sneutrinos may become thermal DM by introducing a new U(l)' 
gauge coupling to make the self-annihilation interaction sufficiently strong [19]. Recently, 
there has been a paper [26] discussing the possibility of right-handed sneutrinos as DM in 
an extended version of the next-to-minimal supersymmetric standard model. 

In the Fp-term hybrid model, a novel possibility opens up. As was first observed 
in [8], there exists a new quartic coupling described by the Lagrangian 2 

4f = l\pN*N*H u H d + H.c. . (3.4) 

This quartic coupling between right-handed sneutrinos and Higgs fields results from the F- 
term of the inflaton field Fs~. \pNiNi + \H u Hd C F$- If strong enough, the interaction (3.4) 
can thermalize the sneutrinos and make them annihilate to a level compatible with the 
current CMB data via the processes depicted in Figure 1. 

For sneutrino masses of our interest, the most relevant processes are the off-resonant 
pair-production of W bosons and the on-shell pair-production of light Higgs bosons. An 
initial estimate of the process NN — > {H u )Hd — > W + W~ for > mw yields 

nDM& .„ (3.5) 
VP A 2 / \ g w m w J 

In order to obtain an acceptable CDM density, relatively large couplings p and A are needed, 
pX > 0.1. However, these large values for A and p are not compatible with the constraints 
derived by inflation. 

The situation differs for sneutrino masses m^ < mw, in large tan /3 scenarios, in which 
light Higgs bosons couple appreciably to 6-quarks [36] . In particular, in the kinematic region 



^or an earlier discussion, see also the paper by N. Arkani-Hamed et al. in [16]. 

2 The implications of a generic singlet-Higgs quartic coupling for the CDM abundance and detection 
were studied before in [27-29], within a simple non-SUSY model. 
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Figure 1: Feynman graphs related to sneutrino annihilation. 

niHi ~ 2m jy , the self-annihilation process JViiVi — >■ {H u )H d — > 66 becomes resonant, and 
the above estimate modifies to 

« lO" 4 x B-\H X - iVxATx) x ( ^g^ ) 2 . (3.6) 

Consequently, if the couplings A, p are not too small, e.g. \p > 1CT 3 , the right-handed 
sneutrino Aq can now efficiently annihilate via a resonant ifi-boson into pairs of 6-quarks, 
thus obtaining a relic DM density compatible with the observed value (1.1). 

We will now show that the naive estimates (3.5) and (3.6) presented in [8] are in a fairly 
good agreement with a complete calculation of all relevant sneutrino annihilation processes 
displayed in Figure 1. To this end, we use the short-hand notation M X y = M(N a N b — > 
XY) to denote the individual matrix elements for the annihilation of sneutrinos N a and 
N b . The contributing processes may be listed as follows (c w = cos0 w , v = 2m w /g w ): 

(i) N a Nb — > H + H~ , via contact quartic interaction and s-channel Higgs exchange: 

3 

M n -, 2 9N a N b H k 9H k H+H~ 

M H+H - = 9 ^ bH+H - - v ^ s _ mk+tmHkTHk ■ (3-7) 
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(ii) N a Nb — ► W + W , via s-channel Higgs exchange: 

2 -,1/2 3 



2 + 1 - 



2m 2 w 



9ff a N b H k 9H k VV 



y — 

s - m H k + %m ^Sn k ' 



(3.8) 



(iii) N a N b — > ZZ, via s-channel Higgs exchange: 

2-i 1/2 3 



M 



ZZ 



g w m w v 

2C,/; 



2+1 



2m| 



9N a N b H k 9H k VV _ 

^ s - + im Hk T Hk ' 



(3.9) 



(iv) iVaiVj, — ► //, via s-channel Higgs exchange: 



M f f = vV2s 

J aja 



with 



3 S/P 
A \- 9N a N b H k 9f a 9 Hk J afa 

s-m 2 Hk +im Hk r Hk 



(3.10) 



(3.11) 



(v) N a N b — > HiHj, via contact quartic interaction, s-channel Higgs exchange and t/u- 
channel sneutrino exchange: 



MHiHj ~ 9N a N b HiHj ~ V ^ ~ m 2 T • TT~ 

2 dNgNcHjdNiNcHj _ 2 9 N a N c H 3 9 N b N c H, 



C=l 



2 ' 
c=l AT C 



(3.12) 



(vi) N a Ni> — > H + W , via s-channel Higgs exchange: 

g w v 



M, 



H+W- 



Am 2 w 



2 i 2 \ 2 

m w + m m 



1/2 3 

9ff a N b H k 9H k H+W- _ 

Z—' s — m?r + ?m#,IV, ' 

fc = l «fc ^fe ^fe 

(3.13) 



(vii) N a Nb — > HiZ, via s-channel Higgs exchange: 



M. 



H t Z 



g w v 
4c, 



.s- / m 2 z + m 2 H \ 2 2 
1 I - 



4m| 



-, 1/2 3 



9N n N h H k 9H k H t Z 



Y — 

fris-m 2 Hk +im Hk r Hk 



(3.14) 
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In the above, the effective sneutrino-to-Higgs couplings gf? a f? b H+H-> 9N a N b HjHj an< ^ 9N a fi c Hi 
that arise from the interaction Lagrangian (3.4) are given by (cp = cos/?, sp = sin/3) 

9N a N b H+H- = -Y c P s f3$ab, (3.15) 
9N a N b HiHj = Y ^ 0(t,ui0 ^ dj + 0ai0 <t>vi s P + OaiO^jCp - ] ai O ] aj S pCp) 

+ (3.16) 

9N a N bHi = y (O^sp + O^iCp) 5 ab , (3.17) 

where O is the 3x3 Higgs-boson mixing matrix, defined such that 

(0 d , U , af = O (#!, H 2 , H 3 f . (3.18) 

For the effective Higgs-boson couplings gu k u,z, gH k H+w-, g^n^^ gH k H+H~, gu k vv and 
gf a g S l l P f 7 , including O, the Higgs-boson masses m# 123 an d their decay widths T Hl23 , 
we follow the notations and conventions of [30, 31] and calculate them by means of the 
computational package CPsuperH. 

The total annihilation cross-section a a b = &(N a Nb — > all) may then be conveniently 
expressed as the sum of all channels, 

3 3 
i=l ij=i f=T,b,t 

The individual cross sections oxy are defined by 

1 1 f t+ 

axY = i i r T7TT7 2 TT / dt\M XY \ 2 , (3.20) 

l + <5 X y 16ttA(s, m% a ,m±J J t - 

with 

1 



x 2s 



TA^^s.m^myA^^.m^m^)) , (3.21) 
X(a,b,c) = (a - b - c) 2 - Abe. (3.22) 

In order to calculate the relic density, we follow [32] and use an effective cross-section 
averaged over all initial sneutrino channels, 

6 

a eS = £ c?ab^(l + A a ) 3 / 2 (1 + A b f/ 2 exp [ - x(A a + A*)] , (3.23) 

a,b=l ^ eS 
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where 

rriM -m<r 



9eS = ]>>(i + A ») 3/ V* Aa , A a = . (3.24) 



In (3.23), both the effects of LSP self-annihilation and co-annihilation with the heavier 
sneutrinos are included 3 . In terms of the effective cross-section (3.23), the thermally- 
averaged effective cross-section may be calculated as 



2-3/2 

(av) 



POO 

/ dvv 2 (a c{i v)e- xv2/4 , (3.25) 
Jo 



2tt 3 /2 

where the integrand is expressed in terms of the relative velocity v, such that 

s = i T77 • ( 3 - 26 ) 

From the expression (3.25), we may determine the freeze-out temperature Xf = m^/Tf by 
iteratively solving the equation 

x f = In TJ ^ TJ - 2 , (3.27) 

V g* Xj / 

where Mpi = 1.22 x 10 19 GeV is the Planck mass and g* is the total number of effective 
relativistic degrees of freedom at the temperature of the LSP freeze-out. The present day 
relic density is then given by 

„ , 9 1.07 x 10 9 Ge^ 1 

^dm h ~ m ( 3 - 28 ) 

J gj mpi 

where J is the post freeze-out annihilation efficiency factor given by 

POO 

J = dvv(<j eff v)etfc(v^x]/2). (3.29) 
Jo 

In our numerical estimates, we neglect the flavor structure of the right-handed sneutrinos 
and treat the three light right-handed sneutrinos iV lj2j 3 as being essentially degenerate 4 . 
Since all three light sneutrinos will contribute to the relic density, we must therefore mul- 
tiply (3.28) by 3 to obtain the final relic DM abundance. 



3 Note that co-annihilation effects become significant, only if the mass differences with the heavier 
sneutrinos are smaller or comparable to the LSP freeze-out temperature, i.e, when — ^ Tf . 

4 Note that the second and third right-handed sneutrinos 7V 2;3 will decay to the LRHS Ni through the 
processes N 2 ,3 ^i7i N\vv. We do not address potential problems for BBN from the late decays of -/V 2j 3, 
since their rates strongly depend on the flavor structure of ptj and the Yukawa couplings [cf. (2.1)] and 
on the details of the model in general. 
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Figure 2: Allowed (mo, 7711/2) parameter space for a mSUGRA scenario with Aq = 300 GeV, 
signp = +, A = p = 10~ 2 and tan/3 = 10 (left panel) and 30 (right panel). The black 
contours show the predicted LRHS mass, while the sneutrino Ni/neutralino Xi/stau fx 
LSP is given by the blue/ green/ orange area. The red area is excluded by direct SUSY mass 
searches. The white contour is defined by the condition rrift — rriH 1 / c 2,, allowing for rapid 
sneutrino annihilation via the H\-boson resonance. 

3.3 Numerical Results 

The numerical analysis is separated in two parts: in the first part, we perform a scan over 
the mSUGRA parameter space to calculate the supersymmetric particle spectrum and 
identify regions where the LRHS can be a possible candidate for CDM. In the second part, 
we specify two mSUGRA scenarios and calculate the constraints on the effective sneutrino 
annihilation coupling Xp by requiring a sneutrino relic density of Odm'* 2 = 0.11. 

In Figure 2 we plot the lightest sneutrino mass as contours in the mSUGRA 
parameter plane (mo^^), for two different values of tan/3 = 10 (left) and 30 (right). In 
both plots of Figure 2, we set A = 300 GeV and p > 0. For the inflaton couplings A, p 
required to calculate the sneutrino masses (3.3), we simply choose 

A = p = 10~ 2 , (3.30) 

in accordance with the bounds (2.26) derived from inflation. The coloured areas in Figure 2 
denote the LSP in the given parameter region: sneutrino Ni (blue), neutralino Xi (green) 
or stau f\ (orange). The red area on the bottom/left is excluded by direct searches for 
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SUSY particles. Specifically, the following experimental mass limits are used [34]: 

104 GeV , 
375 GeV , 

289 GeV , (3.31) 
95 GeV , 
130 GeV . 

Figure 2 was determined by appropriately using universal soft SUSY-breaking parameters 
at the GUT scale according to the mSUGRA scheme and then solving the MSSM RG 
equations down to the electroweak scale. In this respect, our computation was aided by 
the software package SPheno [35] . We neglect the RG running of the sneutrino parameters 
Mj~ and A p which enter the sneutrino mass matrix (3.1), and identify them directly with 
mg and A , respectively. This is a reasonable approximation as their RG evolution is only 
driven by the small couplings A and p. The Higgs coupling parameter p is then calculated 
consistently by requiring proper electroweak symmetry breaking. In the Fo-term model, 
the p term originates from the VEV of the inflaton (2.9). This immediately allows us 
to calculate both the inflaton VEV, (S) = p/\, and the mass scale of the right-handed 
neutrinos, m N = p(S) = (2.1). For the Aq LSP region of interest and with our choice 
A = p = 10~ 2 , p and m N are equal and of order 300 GeV. The mass of the LRHS as 
LSP ranges between 20-100 GeV. This allows for a rapid annihilation of Aq via the Higgs 
resonance, — m Hl /2 57 GeV, along the white contour in Figure 2. 

The Fo-term model puts strong constraints on the mSUGRA parameter space, when 
requiring a sneutrino LSP and taking into account bounds from inflation. As can be 
seen in Figure 2, the connection between LRHS mass Aq and p generally points towards 
a low-energy SUSY spectrum. This coincidentally includes the i?i-boson funnel region, 
where mn x ps ^ m N 1 - On the other hand, very large and small values for A and tan/5 
are disfavoured as they generally exclude a sneutrino LSP. The above correlations may be 
somewhat relaxed if non-universal inflaton couplings A and p are considered. 

In order to compute the sneutrino relic density and analyze the constraints on the 
effective annihilation coupling Ap, the following two mSUGRA scenarios have been selected: 

• Scenario I: 

m Q = 70 GeV, m 1/2 = 243 GeV, A = 300 GeV, tan/3 = 10, p = 303 GeV . (3.32) 

• Scenario II: 

m = 125 GeV, m 1/2 = 212 GeV, A = 300 GeV, tan/3 = 30, p = 263 GeV . (3.33) 



m 



Xi 



> 



rrig > 
m- g > 



vn-n > 



m 9L > 
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Figure 3: Effective annihilation coupling Xp as a function of the mass of the LRHS for 
the observed relic density Qdm^ 2 = 0.11 (blue curve) in the mSUGRA Scenario I (3.32). 
The actual sneutrino and neutralino masses in the scenario are indicated by vertical lines. 
The red curves denote the upper bound on Xp as obtained by the CDMS-II experiment and 
as expected by the projected sensitivities of SuperCDMS and Xenonl T. 



In addition, we keep the LRHS mass as a free parameter. The effective annihilation 
coupling Xp (3.4) is then consistently calculated so as to obtain a sneutrino relic den- 
sity i^DM h? = 0.11, consistent with observation. Furthermore, we assume that the mass 
splitting between the light and heavy right-handed sneutrinos is sufficiently large so that 
co-annihilation can be safely ignored. This is valid as long as there is a sizeable mixing 
between the light and heavy right-handed sneutrino states, which is certainly true for the 
mass range < m^o of our interest. All other MSSM parameters and masses were calcu- 
lated within the mSUGRA framework. Numerical estimates of the allowed parameters in 
the (wjy- , Ap)-plane are shown for Scenarios I and II in Figures 3 and Figure 4, respectively. 

As we have seen in Section 2.3, the requirement for successful inflation puts upper 
limits on the couplings A and p. Given (2.22), (2.23) and 2.26), the upper limits on the 
product Xp for an inflaton sector with a minimal and next-to-minimal Kahler potential 
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Figure 4: As in Figure 3, but for the mSUGRA Scenario II (3.33). 

may easily be deduced to be 

\p < 2.3 x 10" 4 , \p < 5.8 x 10" 4 , 



(3.34) 



respectively, at the soft SUSY-breaking scale M SU sy- On the other hand, Figures 3 and 4 
show that it should be 

\p > 2 x 10" 4 , (3.35) 

in order to account for the observed DM relic abundance in the i^i-boson funnel region, 
where m Larger values of tan/3 do suppress the coupling required to get the 

observed relic density, but not to a level compatible with the inflationary constraints (3.34). 
In general, we find that LRHS masses larger than about 100 GeV are not possible within a 
mSUGRA realization of the F^-term model. This is indicated by the value of the neutralino 
mass in the given mSUGRA scenario as displayed by vertical lines in Figures 3 and 4. 

Further constraints on the (m^ , Xp) -plane may be obtained by taking into account 
the limits from direct searches of experiments which look for scattering between Weakly 
Interacting Massive Particles (WIMPs) and nuclei. Specifically, a WIMP, such as the 
LRHS, can directly be detected through its elastic scattering with a nucleus. In our case, 
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the relevant scattering process is Ni + |J — > iVi + and proceeds via a Higgs-boson 
t-channel exchange. Its cross-section may well be estimated by [29] 

^.nucleus 

(l/2Ap)V|M x | 2 ml d 
* ' { } 

where m re d is the reduced mass of the LRHS-nucleus system, i.e. 

m red = - , (3.37) 

+ m x 

and Mx is the nuclear matrix element. For comparison purposes, we express our results in 
terms of the nucleon cross section. Assuming the nucleus to be composed of A independent 
nucleons, the nuclear cross sections then simply scale quadratically with the nucleon number 
A and the reduced masses: ^r e d(p) a ci uclcus = ^ 2m red(z^)°"ei ucleon - The nucleon matrix 
element M nuc i eon ~ 10~ 3 is mostly sensitive to the strange-quark Yukawa coupling. An 
adequate estimate of the elastic scattering cross section <j° ucleon of a right-handed sneutrino 
with a nucleon yields [29] 



nucleon ^ ^ m-50 
cl 



Xp V f 10Q GeV\ 4 /50 GeV 



(5xl0-o^)(^ (^1) . (3,8) 



m N 1 



The upper limits on Xp are derived by comparing the estimate (3.38) with the current 
bound on the spin-independent nucleon cross section from the CDM S-II experiment and the 
expected sensitivities of the SuperCDMS extension [37] and the XenonlT experiment [38]. 
These limits are included in Figures 3 and 4. The current bound already excludes large 
parts of the (rrift , Xp) -parameter plane, except of the Higgs-boson funnel regions. In the 
near future, the upgraded experiment SuperCDMS will cover a large part of the parameter 
space, but it will leave open the lightest Higgs-boson pole region which is theoretically 
favoured by inflation within the mSUGRA framework. The proposed XenonlT experiment 
is expected to further narrow down this uncovered parameter range of the Fo-term model. 

Dark Matter may also be indirectly searched for through the detection of its final 
annihilation products, such as photons, positrons, anti-protons or neutrinos. The domi- 
nant channel of the LRHS annihilation in the Higgs funnel is determined by an effective 
scalar coupling with a bb pair, which is approximately independent of the relative veloc- 
ity of the annihiliating sneutrinos. Rates at low temperatures resulting in gamma-ray or 
charged particle fluxes are therefore not suppressed compared to the rates at the freeze-out 
temperature responsible for the LRHS relic density. There are several signals that could be 
explained as an observation of DM annihilation but, as of now, do not provide a consistent 
picture interpretable by a single DM candidate and model. For example, the excess in the 
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diffuse galactic gamma ray spectrum measured by the EGRET detector may be interpreted 
by a 50-100 GeV WIMP, as given by the LRHS in our model, whereas the 511 keV line 
observed by the INTEGRAL satellite would hint at an MeV DM particle [39,40]. Up- 
coming projects such as the GLAST and PAMELA satellites will have higher sensitivities, 
probe new energy ranges and should provide a clarification of the observational status. 
High-energy neutrinos as annihilation products are expected and can be searched for in the 
Sun and the Earth, as WIMPs can accumulate in their centre. For the LHRS there is no 
spin-dependent coupling to nuclei, and its capture rate along with the produced neutrino 
flux is suppressed. In addition, for an annihilation via the Higgs resonance, the effective 
annihilation coupling required to get the correct relic density is very small. The LRHS is 
therefore not expected to be within the reach of high-energy neutrino telescopes [41], such 
as IceCube [42]. 



4 Conclusions 

We have analyzed in detail the relic abundance of the lightest right-handed sneutrinos 
(LRHS) in the supersymmetric F^-term model of hybrid inflation. The inflationary po- 
tential of the model results from the F-term of the inflaton multiplet S. The F^-term 
model also includes a subdominant non-anomalous D-term generated from the local U(l)x 
symmetry of the waterfall sector, which does not affect the inflaton dynamics. As was 
mentioned in the introduction and further discussed in Section 2, the model adequately fits 
the current CMB data of inflation and provides a natural solution to the so-called gravitino 
overabundance problem, without resorting to an excessive suppression of possible renor- 
malizable couplings of the inflaton to the MSSM particles. Finally, the F^-term model 
closely relates the /x-parameter of the MSSM to an SO (3) symmetric Majorana mass m N 
through the VEV of the inflaton field. If A ~ p, this implies that /i ~ m N) so the model may 
naturally predict lepton-number violation at the electroweak scale and potentially account 
for the BAU via thermal resonant leptogenesis. 

In spite of the explicit lepton-number violation through the Majorana term | p SNiNi, 
the Fo-term hybrid model conserves F-parity. Consequently, the LSP of the spectrum is 
stable and so qualifies as candidate to address the CDM problem. The new aspect of the 
Fo-term hybrid model is that thermal right-handed sneutrinos emerge as new candidates 
to solve this problem, by virtue of the quartic coupling: \ \p N*N*H u H d + H.c. This new 
quartic coupling results in the Higgs potential from the F-terms of the inflaton field, and 
it is not present in the more often-discussed extension of the MSSM, where right-handed 
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neutrino superfields have bare Majorana masses. Provided that the couplings A and p are 
not too small, e.g. A, p £ 10~ 2 , the LRHS A^sp as LSP can efficiently annihilate via the 
lightest Higgs-boson resonance Hi into pairs of 5-quarks, in the kinematic region mu 1 ~ 
2m^ Lgp , and so drastically reduce its relic density to the observed value: fi DM h 2 0.11. 

Experiments, such as CDMS-II, SuperCDMS and XenonlT, which look for signatures 
of WIMPs through their elastic scattering with nuclei, will significantly constrain the al- 
lowed parameter space of the Fo-term model. They will exclude most of the parameter 
space, except possibly of a narrow region close to the lightest Fi-boson resonance, where 
mux ~ 2 m N L p - K might seem that to obtain this particular relation between the masses of 
the Hi boson and Alsp, a severe tuning of the model parameters is required. However, it is 
worth stressing here that such a mass relation may easily be achieved within a mSUGRA 
framework of the F^-term model that successfully realizes hybrid inflation. 

The LRHS scenario of CDM requires relatively large A and p couplings that could, 
in principle, make Higgs bosons decay invisibly, e.g. H — > Alsp A"lsp- Also, right-handed 
sneutrinos could be present in the cascade decays of the heavier supersymmetric parti- 
cles. The collider phenomenology of such a CDM scenario lies beyond the scope of the 
present article. Instead, we note that the F^-term hybrid inflationary model can give 
rise to rich phenomenology which can be probed at high-energy colliders [20,21], as well 
as in low-energy experiments of lepton flavour and number violation, such as Qv(3(3 de- 
cay, p — > e7 [43], p — > eee and p — > e conversion in nuclei [44,45]. It would therefore 
be very interesting to systematically analyze possible correlations between predictions for 
cosmological and phenomeno logical observables in the Fo-terni model. 
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